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Abstract 

3 We analyze the motion of individual beads of a polymer chain using a 

discrete version of De Gennes' reptation model that describes the motion of a 
polymer through an ordered lattice of obstacles. The motion within the tube 
can be evaluated rigorously, tube renewal is taken into account in an approx- 
imation motivated by random walk theory. We find microstructure effects to 
be present for remarkably large times and long chains, affecting essentially 
i . all present day computer experiments. The various asymptotic power laws, 

C**^ | commonly considered as typical for reptation, hold only for extremely long 

chains. Furthermore, for an arbitrary segment even in a very long chain, we 
t— I ■ find a rich variety of fairly broad crossovers, which for practicably accessible 

chain lengths overlap and smear out the asymptotic power laws. Our anal- 
. ysis suggests observables specifically adapted to distinguish reptation from 

motions dominated by disorder of the environment. 
£3 key words: reptation, polymer dynamics 
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O ■ I. INTRODUCTION 

O ■ 

. £h ! The reptation concept, introduced long ago by De Gennes |1| , has become a basic scenario 

of polymer dynamics. It deals with the motion of a flexible chain surrounded by a network 
of impenetrable obstacles. Originally it has been proposed for understanding the diffusion 
of a polymer through a gel, but it widely is used also to analyze the dynamics of polymers 
in a melt or in a dense solution. A recent review may be found in ||. 

The basic idea is best explained for polymer diffusion through a gel that is swollen by 
some solvent. The gel consists of polymer strands that are chemically linked together to form 
a three dimensional network. In the solvent filled space between the strands some free chain 
not chemically linked to the gel may exert a Brownian motion. The hard core interaction 
with the strands of the gel constrains this motion to a 'tube' that is defined by the strands 
surrounding the instantaneous chain configuration. Due to crosslinking, these 'topological 
constraints' are essentially fixed. (The thermal motion of the gel is neglected.) The array of 
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topological constraints forming a tube around the instantaneous chain configuration is one 
important ingredient of reptation. The second important ingredient of De Gennes' theory 
is the specific mechanism of microscopic chain motion within the tube. Since the chain is 
highly flexible, it cannot move as a whole like a rigid body, but it moves through thermal 
displacements of individual segments. It also does not lie stretched in its tube, but as a 
consequence of flexibility and entropy, there will be a finite density of small wiggles or side 
loops of 'spared length'. These wiggles, called 'defects' by De Gennes, exert a diffusive 
motion along the chain and thus along the tube. This is taken as the basic dynamical 
mechanism. 

If a defect reaches a chain end, it can be destroyed, thus prolonging the tube by its spared 
length in some random direction. A chain end can also contract into the tube, thus creating 
a defect and destroying the end of the previous tube over a distance of its spared length. 
This mechanism results in a thermal motion of the chain ends, which gradually destroys the 
original tube until finally the chain has found a completely new configuration. 

Accordingly, the calculation of dynamical correlation functions in the reptation model 
can be decomposed into two well separated steps. The first step is the one-dimensional 
diffusion of conserved defects along the chain as the basic thermal motion. This yields the 
curvilinear displacement of some interior segment along the tube, and tube renewal at the 
chain ends. The second step consists in embedding the segment displacement into real space 
by taking the (equilibrium) random walk configuration of the chain and thus of the tube 
into account. 

Treating defect dynamics in a continuum approximation, De Gennes U showed the 
existence of several characteristic time scales and derived asymptotic power laws. We here 
recall some of the results, using notation adequate for a discrete chain composed of N 
segments. 

There emerge three characteristic time scales. A microscopic time 



is needed for the chain motion to feel the constraining environment. For shorter time inter- 
vals the chain does not properly realize that it is confined to a tube. A second time 



is of the order of the Rouse time, i.e., the longest relaxation time of a free chain. In the 
present context it, however, has a quite different interpretation: it is the time a defect needs 
to diffuse over the whole length of the chain. The third scale 



is known as 'reptation time'. Time intervals of the order of T 3 are needed to completely 
destroy the original tube. 

We now consider the average displacement in space of some internal bead j of a chain of 
length A^: 



T 



o 



T 2 ~ N 2 



T 3 ~ N 3 




(1.1) 
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The bar indicates the average over all defect motions along the chain, and the pointed 
brackets denote the average over all tube configurations. The reptation model predicts 



gi(j,N,t) 



{t/Nf/ 2 T 2 < t < T 3 
t/iV 2 T 3 <t . 



(1.2) 



Note the asymptotic nature of these results. In precise mathematical terms the first line of 
Eq. (1.2), for instance, reads (we use the notation j = j/N): 



These predictions therefore concern long chains and times large on microscopic (T ) or even 
mesoscopic (T2,T 3 ) scales. Similar results for other quantities like the motion of the center 
of mass have also been derived. 

Much work has aimed at extracting these power laws from simulation data or physical 
experiments on polymer diffusion mostly through melts, dense solutions, gels or spatially 
fixed, but disordered configurations of obstacles. On a superficial level, the expected behavior 
is seen, but closer inspection proves the outcome to be inconclusive. For instance, for the 
central bead of the chain one finds some slowing down of the motion that according to most 
simulations in some intermediate range seems to follow an effective power law 



This regime immediately is followed by a smooth crossover towards free diffusion: g\ ~ t, 
with no proper indication of a t 1//2 -regime. For the center of mass motion the results likewise 
are unsatisfactory. 

The interpretation of these results is quite ambiguous. In melts and dense solutions a slow 
relaxation of the surrounding medium, known as 'constraint release', will play a role ||. For 
gels, there exist hints that some types of disorder may strongly modify the chain dynamics 
^-0]. Also a quite trivial effect should be noted: The bead motion of a simple Rouse chain in 
a straight tube of finite diameter shows a crossover from three dimensional motion at short 
times to one dimensional motion at larger times, that closely resembles simulational results 
found in the literature, but that has little to do with reptation. Nevertheless two results of 
the simulations deserve special interest here: Following the change of chain configurations 
with time, Kremer and Grest |J convincingly demonstrated that the tube concept is valid 
even in melts. Again in a simulation of melts, Shaffer || did a controlled study of the 
influence of chain impenetrability, showing that it strongly modifies the motion of a test 
chain. 

This situation clearly longs for a detailed study of the pure reptation scenario, suppressing 
all additional background effects. A first notable attempt was undertaken by Evans and 
Edwards |TIJ. Their work, however, suffered from short chains and poor statistics. In their 
model they considered the motion of a random walk chain on a cubic lattice with a second 
(dual) cubic lattice (with lattice points on the centers of the cubes of the original lattice) 



lim 



1//4 lim gi(j,N;t) = const 



j fixed 



OG 



N — > oo 
t, j fixed 
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representing the obstacles. The two-dimensional version of the model is shown in Fig. 1. We 
have carried through new extensive simulations of this model, and our results are described 
in detail in a second paper III]. Here we only note that we indeed find a t 1//4 -regime, (as 
well as indications of a t 1//2 -regime) , but even under the most favourable circumstances it is 
fully developed only for quite long chains and surprisingly long times ]I2] . Most of the time 
regime is covered by broad crossover regions. 

To understand these results, it needs a more detailed analytical evaluation of the repta- 
tion model, giving quantitative results for functions like gi(j,N,t) also outside the asymp- 
totic limits. Particular attention should be paid to the time regime t < T 2 , since present day 
simulations for long chains hardly can reach the regime t > T 2 . Furthermore we should look 
for quantities which allow for a clear distinction between reptation and Rouse-type motion 
in a disordered environment. The present work is devoted to that task. 

In view of the (computer-) experimental situation we also should account for short time 
effects. Though these effects to some extent will depend on the microscopic model, it for the 
interpretation of the data is most important to get some feeling for their typical character 
and range. We thus use a fully discretized version of the original De Gennes model, adequate 
for an ordered lattice of obstacles. Depending on the time regime we calculate g\ (j, N, t) and 
related quantities either rigorously or in a seemingly quite good approximation. For large 
times some of our results with small modifications also could be extracted from equations 
given in Jl|. As shown in the second paper [|ll|] our analytical results perform very well in 
explaining the Monte Carlo data. They thus provide a firm basis for the study of polymer 
motion through more complicated environments like melts or a disordered background. A 
short account of some of our analytical and numerical results has been published in ||12|| . 

We should delimit our work from two related, but different lines of approach. Rubinstein 
|T3| proposed a similar lattice model, that nowadays widely is used in the theory of elec- 



trophoresis. Compared to his model our approach is simpler, since without external fields the 
defect motion decouples from the chain configuration. As a result, we can calculate dynamic 
correlation functions analytically. Doi and Edwards fll4| , in their work on melts, retained 
the tube concept but abstracted from the specific dynamical mechanism of defect motion. 
They consider the diffusion of a so called 'primitive chain', which is a fictitious object fixing 
the tube configuration. In the large time regime t > T 2 , this approach is equivalent to 
De Gennes' model, but for shorter times the intrinsic dynamical mechanism should not be 
ignored. Doi |15| extented this theory by adding internal Rouse dynamics, thus considering 
a (stretched) one dimensional Rouse chain in a tube. If we ignore the discreteness of the 
defect motion, our results reduce to those of this model, and for those aspects of the tube 
renewal specifically considered in |T5[, we then find similar results. 

This paper is organized as follows. In Sect. II we define our model, and we derive expres- 
sions for physical observables. We mainly will consider the single-bead motion. Quantities 
like the center of mass, which involve two-bead correlations, in principle are accessible also, 
but will not be considered here. In Sect. Ill we derive the generating function of the mo- 
ments of the single bead motion. It yields rigorous expressions for quantities like g±(j, N, t), 
as long as the bead j on average still sits in a part of the original tube. The analytical and 
numerical evaluation of our results in that time regime is contained in Sect. IV. Sect. V deals 
with end effects, presenting a step towards the description of tube renewal. The calculation 
of gi(0, N, t) amounts to determine the average maximal value within the time interval [0, t] 
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of a correlated stochastic process. For this problem no rigorous solution is known. We use a 
'mean hopping rate' approximation, that is rigorously valid for uncorrelated, i.e., Markovian 
processes. In the same spirit, we then evaluate the effect of tube renewal on some arbitrary 
bead j for times t < T 3 , and we estimate the reptation time T 3 . Finally Sect. VI summa- 
rizes our findings. Some useful formulae and more complicated calculations are presented in 
appendices. 



II. MODEL AND OBSERVABLES 
A. Construction 

In constructing our model, we are guided by reptative motion of a lattice chain. Assume 
for simplicity that the chain is represented by a random walk of N steps (segments) on 
a square lattice of lattice constant £q. The configuration is fixed by the bead coordinates 
{r (£), • • • , r N (t)}, which give the positions of the N + 1 endpoints of all the segments. A 
lattice of impenetrable obstacles occupies the centers of the cells of the first lattice, see Fig. 
1. The obstacles exclude any kink jumps or more complicated modes of motion, and leave 
only tips of hairpins, i.e., beads j with r-j +1 = r,-_i, free to move, Tj then can move between 
all sites neighboring r,-_i. If the hairpin happens to lie down on the chain, it can stand up on 
a different bead, and so diffuse along the chain. If a hairpin (i.e., a 'defect') diffuses across 
a bead, it displaces it by the spared length 2£ along the chain configuration. Of course, 
also the endpoints are free to move, and at the ends hairpins are created or destroyed. In 
our simulations we used the three-dimensional version of this model, that originally was 



proposed by Evans and Edwards JlQ 



Abstracting from this concrete realization, we consider the chain configuration as a ran- 
dom walk of N segments of fixed length that henceforth will be taken as the unit of 
length. The defects are modelled as particles sitting on the beads of the chain. Each parti- 
cle represents a defect of spared length i s . It can hop to each of the two neighbouring beads 
along the chain with probability p. For a given particle, the other particles are assumed to 
be just a part of the chain, so the particles do not interact, and each particle diffuses not 
along a primitive chain, but along the whole random walk configuration of the chain with 
length N, which corresponds to a physical chain of length N + £ s . This allows us to take 
independent averages over the particle diffusion along the chain and over the random walk 
configuration of the chain in embedding space. We introduce the equilibrium density po of 
particles as the average density of particles per bead. The assumptions of the model will be 
further discussed in Subsect. II. B below. 

As pointed out above, at the chain ends (beads and N) defects can be destroyed or 
created. We can turn this grandcanonical problem into an easier treatable canonical one 
by coupling the chain ends to large virtual particle reservoirs, located at bead or N, 
respectively. A particle may hop from bead 1 (N — 1) into reservoir (A) with probability 
p. It may hop out of the reservoir onto the chain (from to 1, or from A to A — 1, resp.) with 
probability p res . We put M particles into the system. Taking p res — > and M — > oo, with 
ppo = PresM/2 fixed, each reservoir becomes a source of particles hopping independently 
into the chain and also a sink absorbing particles from the chain end, such that the proper 
particle dynamics at the chain ends is achieved. This canonical formalism allows a very 
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simple calculation of the displacement of some bead through simple book keeping of the 
change of particle number on one side of the bead. 

The calculation of the displacement of a bead j thus decomposes into three steps: (i) 
Calculate the motion of an individual particle on the one-dimensional chain, (ii) Calculate 
the balance of particles, that have moved over the bead j within the time of considera- 
tion. (Hi) Determine from this balance the displacement of the bead in three-dimensional 
embedding space. 

Let us first characterize the motion of the individual particles on the chain. Pi(j,t) is 
the probability, that a given particle at time t sits on bead j, < j < N. This probability 
develops in (discretized) time according to the equation 

N 

Pi(j,* + 1) = E W n' P iU',t) , (2.1) 

j'=0 

with the hopping matrix for all < j < N 

W jr = (1 - 2p)6 jj , + p(5 jJ/+1 + S jd ^) for 1 < f < N - 1 

W j0 = Pres^jl + (1 - Pres)8j0 

W jN = Pres^jN-l + (1 - Pres)5jN ■ (2.2) 

The normalized stationary solution of Eq. (2.1) takes the form 

P[ eq \j) = {l/p)/[{N-l)/p + 2/p res ] 

A N (0) = P[ eq \N) = pPl eq \j)/p res (2.3) 

for 1 < j < N — 1. We have defined the equilibrium density of particles on the chain (beads 
1 to N — 1) to be p , and we thus have for any one of the M particles 

P[ eq \ 3 ) = ^, 1< 3 <N-1. (2.4) 

This fixes p res 



rn , — " n i — i 

and results in 

P (e,) (0) = p m {n) = i ^ _ n__± po + Q_yj (2 6) 

So far we have defined a standard and readily solved problem of hopping dynamics. The 
special features of reptation emerge when we turn to the physical observables, filling in steps 
(ii) and (Hi). We first consider the motion of a bead j within a time interval short enough 
that it stays in the tube defined by the initial configuration {rj(0)}. (See Fig. 2.) At time t 
the displacement of that bead along the tube is given by 

e a n(j,t)=e 8 [n-(J,t)-n + (J,t)] , (2.7) 
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where n± (j, t) is the number of particles that have passed over bead j within time interval 
[0,t] from the 'left' (j' < j) or from the 'right' (j' > j), respectively. (Thus n(j, 0) = 
by definition.) Since the tube configuration itself is a random walk, bead j in real space is 
displaced by a random walk of £ s \n(j,t)\ steps. We thus find 



g 1 (j,N,t) = ((T j (t)-T j (0)y 



= £ s \n(j,t)\ . (2.8) 

We here and in the sequel pay no special attention to the situation where at time or t a 
particle sits right on top of bead j. In comparison to experiments such configurations will 
yield a small time independent contribution to g±(j, N, t), (which we discuss further in [0). 

We now turn to the motion of the ends of the chain. This intrinsically is the problem 
of tube renewal. We again consider the situation of Fig. 2, where at time t a part of the 
original tube still exists that at time was formed by the beads j< < j' < j > . (Note that j 
is an index of the bead at all times, and of the tube at time t = 0.) We concentrate on the 
motion of chain end j = 0. The length j < on which the original tube is destroyed from the 
left till time t, is determined by the minimal particle occupation of reservoir j = within 
time interval [0,t], or equivalently by the maximal contraction of the chain into the tube 

j< = L n max (t) , (2.9) 



n max (t) = max {n(0, s)} > . (2.10) 

se[0,t] 

The newly formed piece of the tube at time t has length j< — £ s n(0, t). Thus the end bead 
effectively has been displaced by a random walk of length 2j< — l s n(0, t), and we find 

^(0,iV,t) = ((ro(t)-r (0)) 2 ) 

= 2£ s n max {t)-£ s n{0,t) . (2.11) 

Since the particle motion obeys detailed balance, the stochastic process n(j, s) is invariant 
under a change of sign, and the last contribution in Eq. (2.11) averages to zero. Our final 
expression reads 

gi (0,N,t) = 2£ a n max (t) . (2.12) 

An end bead is immediately subject to tube renewal. If we now consider some internal 
bead j, < j < y, it initially will move within the original tube, but will feel the tube 
renewal for sufficiently large times. In discussing this superposition of the two modes of mo- 
tion we assume t < T 3 , such that a part [j<, j>] of the original tube still exists. Furthermore 
we restrict ourselves to j y, so that bead j is affected by tube renewal from the closer 
chain end only. For a given realization of the stochastic process we have to distinguish 
two situations, corresponding to the two cases discussed above: 

a) j< < max{ j,j + £ s n(j, £)}, 

i.e., either the initial position or the end position of bead j at time t is still in the 
original tube. Then the bead is displaced by a random walk of £ s \n(j,t)\ steps. 
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b) j< > max{j, j + £ s n(j, t)}. 

The bead is displaced by a random walk of (j < — j) + (j< — j — £ s n(j, £)) steps. 

Combining these results and taking averages we find (recall Eq. 2.9) 



9i(J,N,t)=e a \n(j,t)\ [l-e(£ s n max (t)-j)e(£ s n max (t)-j-£ s n(j,tj)] (2.13) 
+ 4 (2n max (t) - n(j, t) - j-^j ®(t s n max {t) - j) Q{i s n max (t) - j - £ s n(j, t)) , 

where we introduced the discrete 9-function 

It is easily checked that this expression generalizes the special cases discussed above. Eqs. 
(2.8), (2.12), (2.13) are exact expressions in our reptation model. They hold as long as a 
part of the initial tube survives, i.e., for times t < T 3 . 



B. Discussion 



Our model is no one-to-one transcription of any model used in simulations. Some modi- 
fications have been mentioned above: The particle moves along an effective chain of length 
N, the physical chain being longer by the spared length of the corresponding defect. If a 
particle sits right on top of bead j, the spatial position of the corresponding bead of the 
physical chain may differ from rj by a vector of length l s j1. Also end beads can move around 
without forming a defect, a motion not taken into account. All these effects, however, are 
relevant for the motion on microscopic time scales only, and they can easily be corrected 
for in comparison with experiments. Furthermore our model implicitly assumes a tube of 
microscopic diameter ~ £ , leading to pT I. We will see, however, that the discrete 
nature of the motion is felt for much longer times. 

The assumption of noninteracting particles may be more serious. Even in a random 
walk chain, realistic defects experience a kind of interaction, if they meet on the chain. 
They may merge and form a larger defect, which has to decay before the defects can travel 
on. This corresponds to a scattering process in the particle picture. In our model the 
average effect of these processes is absorbed into the effective hopping rate, which is a model 
parameter that sets the elementary time scale. Only anomalously large fluctuations in the 
local defect density could invalidate our results. Since we are concerned with the dynamics 
under equilibrium conditions, we can estimate the probability of such fluctuations from the 
equilibrium configuration of a random walk on a hypercubic lattice. For a simple estimate 
note that the density of defects is proportional to 1/ (2c?) with 2d the coordination number 
of a hypercubic lattice. The probability of larger defect structures is found to decrease 
exponentially with the number of defects involved. We therefore expect that the scattering 
processes do not seriously affect the model. 

We also would like to comment on two questions related to the basic philosophy of the 
approach. Firstly, where in our model does the assumption of a regular lattice of obstacles 
come into play? This happens in two ways. If the obstacles are irregularly distributed in 
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such a way that they occupy some spatial volume otherwise accessible to the polymer, the 
polymer or parts of it might be confined to entropic traps. As a consequence the simple 
random walk type embedding of the tube into real space is invalid. A second effect of disorder 
is present even for purely topological constraints, which hinder the motion but exclude no 
finite space for the static polymer configuration: In such an environment, the random walk 
configuration of the chain is unchanged, but the mobility of the defects depends on their 
spatial position. Both effects effectively couple the defect motion to the chain configuration. 

Secondly, we comment on our notion of the tube. In a lattice model like that described 
at the beginning of this section and in Fig. 1, the tube may be defined as the non- reversal 
random walk configuration of the chain constructed by cutting off all hairpins, double hair- 
pins, etc. Since, however, defects have to move along the whole chain, not only along its 
backbone, the thus constructed tube has no immediate relevance for the microscopic mo- 
tion. We thus here use the word 'tube' as synonymous to the instantaneous full random 
walk configuration of the chain. This definition allows us to proceed much further in our 
analytical approach. 

We finally comment on the generalization of our model to self-repelling (excluded volume) 
chains. First of all, if the reptation model is applied to melts or dense solutions, it is typically 
argued, that the self-repulsion of the individual chain is screened by the presence of the other 
chains, such that each chain, though repelling, has a random walk configuration over longer 
distances. In an regular array of fixed obstacles, however, self-repulsion has nontrivial effects. 
Technically, we have to distinguish the effect of the repulsion among beads spaced a small 
or a large distance along the chain. On small distances, \ji — J2I = 0(1), the repulsion 
among beads j±, J2 will change the precise physical realization of the defects as well as their 
density, and it will lead to a readjustment of the effective hopping rate, similar to the 
scattering processes addressed above. This consideration also applies to models for melts 
etc. The nontrivial aspects of the excluded volume problem, however, are due to interactions 
of beads spaced far along the chain. Now it is known that the excluded volume interaction 
decreases the probability of close encounter of such pairs of beads, to the extent that the 
number of close encounters, appropriately corrected for the effects of pairs close along the 
chain, vanishes for an infinitely long chain [|T(J . The excluded volume therefore should have 
no nontrivial influence on the defect dynamics. Only the static configuration of the chain 
must respect excluded volume statistics. We thus conclude that Eq. (2.7), for instance, 
for an excluded volume chain in a regular obstacle lattice to first approximation should be 
replaced by 

9i(j,N,t)=£ 3 (\n(j,t)\) , (2.15) 

where v m 0.588 is the excluded volume (Flory) exponent. A more accurate treatment 
should build upon the known internal correlation functions of the excluded volume problem. 
Our other expressions are subject to similar modifications. 

III. GENERATING FUNCTION AND MOMENTS 

To calculate the thermal displacement gi{j,N,t) in space of a bead j of the chain, we 
have to control the stochastic variables n(j, t) and n(0, s) for 1 < s < t of particle motion 
along the chain. We thus introduce the distribution function 
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f > ({n s },n;j,t) = 5 nn{j:t) ]]_ 5 nan{0 ,s) , (3.1) 



s=l 



as well as the generating function 



V({x s },x;j,t) = ]T Y, x n Hx^V({n s },n;j,t) 

n=-oo {n s }=-oo s=l 



t 

= x n ^t) Y[ x n (°' s ) . (3.2) 

s=l 

To find a manageable expression for n(j,t), we now take advantage of the fact that by 
construction of the model, the number of particles is conserved. We thus can calculate n(j, t) 
as the change in the number of particles found at positions f < j: 



M 

n(j,t)= ]T [e(j-j(m,0))-e(j-j(m,t))] . (3.3) 

m=l 



This expression is valid for all j including j = and N. j(m, s) is the bead where the m th 
particle sits at time s. Eq. (3.3) thus expresses n(j,t) in terms of single-particle coordinates. 
Since the M particles move independently, we now can rewrite Eq. (3.2) as 

V({x s }, x- j, t) = {Q({x s }, x- j, t)] M , (3.4) 

where Q{. . .) is a single particle contribution. 



Q({x s }, x- j, t) = x®ti-io)-eU-h) 7J x *>*> (3.5) 



s=l 



(Note that from Eq. (2.14), O(0 — j s ) = 5oj s , j s > 0). Here j s is the position on the chain 
where the particle considered sits at time s. Writing out the average explicitly, we have to 
evaluate the expression 



N t 



Q({x s }, x; j, t)= J2 xW-M-W-*) H ^o-^ s 

{jo,ji,-,jt}=0 s=l 

■ A Pi eq) Uo) (3.6) 

s'=l 

Aiming at the limit M — > oo, we need to calculate Q(. . .) including terms of order 1/M. 
Since both the hopping rate p res from a reservoir onto the chain and the probability P\ eq \jo) 
with 1 < jo ^ N — 1 to find the particle on the chain initially, are of order 1/M (cf. Eqs. 
(2.4), (2.5)), we need to consider two types of processes only: 

i) The particle sits on the chain at time t — 0. If during its walk it is caught by a 
reservoir, it is not allowed to leave it again. 

ii) At t — the particle sits in a reservoir. It may leave it, but if its caught by a reservoir 
again, it may not jump on the chain a second time. 
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The contributions of these processes are easily expressed in terms of the Green's function 
for the particle walking on the chain. 

G$(8) = (W'). r , l<j,j'<N-l (3.7) 

Here Wj? is the matrix W (Eq. (2.2)), restricted to the subspace 1 < j,j' < N — 1. To 
give an example we consider a special process of type ii): namely that the particle sits in 
reservoir at time 0, and on some bead f with j < f < N at time t, i.e., still on the chain, 
but on the other side of bead j. The contribution of all such processes to Eq. (3.5) reads 

t N-l t 

E E x ft x^P^WiWoo) 3 - 1 Pres Gf l (t-s) = 

s=1 j'=j+l s'=s 

s=l j'=j + 1 s'=s 

where in the last line we used Eqs. (2.2), (2.5), (2.6). The factor of x accounts for the fact 
that jo — < j, whereas jt = j' > j. Furthermore, after the particle has left the reservoir, 
each time step s' contributes a factor x s >. All other possible processes can be evaluated 
in the same way. The calculation is straightforward. The expressions can be simplified by 
using the well known explicit form of G^y(s): 

^) = ^ K f™(fj)™(f ?)°t (3-8) 
a fe = l-4psin 2 l— j (3.9) 



We here only summarize the final results for (3.5): 

1 

M 



Q({x s },x;j,t) = l + -]-\Q 1 (x-,j,t) + Q 2 ({x s y,t) + Q 3 ({x s },x-,j,t)} , (3.10) 



Q 1 (x;j,t) = p (x+ l - - 2^ A 1 (j,t) , (3.11) 
Mj,t)- N + 2N ^(1 a k ) s . n2 ^ , (3.12) 

Q2({xs};t) = p oP j2[i[ xj+fl ^'-2j (3.13) 

s=l \s'=s s'=s / 

+pop 2 e ( n ^ s - 1 ) ( n x; 1 -!) a 2 ( S2 - Si ) , 

0<s 1 <s 2 <t \s=si I \S=S2+1 I 



Ms) = j^Y, «^in 2 ^ ) , (3.14) 



k=i 
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X,' 




I irk 



N N 



cos 



fc=l 



N 



] + 



(3.15) 
(3.16) 



Some useful properties of the coefficients A b , b = 1,2, 3, are collected in appendix A. Substi- 
tuting Eq. (3.10) into Eq. (3.4) and taking M — > oo we find 



V({x s },x;j,t) = esxp{Qi(x;j,t) + Q 2 ({x s }; t) + Q 3 ({x s },x; j,t)} 



(3.17) 



We have written this result in a form that factorizes into a contribution Q\ giving the 
statistics of n(j,t), a contribution Q 2 giving the statistics of the occupation {n(0, s)} of 
reservoir 0, and an interference term Q 3 . Note that summation over some variable n s 
or n in V({n s },n; j,t), simply amounts to setting the corresponding variable x s or x in 
V({x s },x;j,t) to 1. 

At this point we can formally establish the connection to the model of a Rouse chain in 
a tube fT5[ . We write x = e lip , x s = e tlfs and we expand to second order in <p, cp a to find 

Ql + Q2 + Q3 Q<p({<Pa},<P'J,t) = po V 2 Sjj{t) 

- y E V^iV 9 ^ (£oo(si) + ^00(^2) - ^oo(|s 2 - si|)) 

si,s 2 =l 

t 

- pop ¥s (Goj(s) + Goj(t) - Goj(t - s)) , (3.18) 



8=1 



where 



Gjf(s) 



ps 1 
N + 2N 



N-l 



E(i- 



k=l 



COS 






) cos 








sin 


2 f irk \ 
\2N) 





(3.19) 



We here used Eqs. (A. 7), (A. 12) of appendix A to evaluate summations over A 2 (s) and 
A 3 (j, s). In the limit s — > 00 and fixed, Gjf(s) reduces to the propagator of the bead 
motion in a one dimensional Rouse chain of length N — 1. For the moments which can be 
derived from Q^, we in this limit recover the results of the model of a Rouse chain in a tube. 

Returning to our full expressions we note that the simplest quantity of interest to us 
in the sequel, is the distribution of n(j,t), irrespective of the occupation of the reservoirs. 
Putting x s — 1 for all 1 < s < t, we find 

1 f dx 

Vi(n;j,t) = — f — x~ n exp(Q 1 {x;j,t)) . (3.20) 
2m J x 

Qi(x;j,t) (Eq. 3.11) is invariant under the substitution x — > 1/x, which leads to 

V 1 (n;j,t)=V 1 (-n;j,t) . (3.21) 
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This expresses detailed balance to hold in our equilibrium problem. (Of course, also invari- 
ance under the substitution j — > N — j holds, which reflects the symmetry of the chain.) 
With the substitution x = e tip , we may express Vi(n;j,t) in terms of a modified Bessel 
function of first kind I n (z): 

VMht) = e- 2 ^ A ^ 7 n (2p A0',0) • (3-22) 
For large values of poA^j, t), we find n 2 ~ p ^i(j, t), and V\ tends to a Gaussian 

A(n; - (47rp A 1 (j,t))- 1 /2 exp (^-—^—^ . (3.23) 

We observe that this limit is equivalent to the quadratic approximation (3.18) discussed 
above. 

With the explicit form (3.11) of Qi, we immediately can calculate all the cumulants of 
our distribution 



d_\ 
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(n(j,t)) 2£ = \x —j In Vi(x;j,t) = 2p A 1 (j,t) 

= ^0*). (3.24) 



The cumulants (n(j, t)) 2e are thus independent of £ for all integer t. Odd moments, of course, 
vanish. We, however, are interested also in odd powers of absolut values. Specifically we 
need to calculate 



+oo 



\ n 0,t)\ = M *Pi(n;j,t) 

n=— oo 



1 / dx eXp[gi(x ;f )] (3.25) 
iri J\x\<i {x - I) 2 



With the substitution x — | + 1 — \ \Jy{y + A) it is a standard exercise to reduce this 
expression to the form 

Rj^=^= (p A 1 (j : t)) 1/2 [l-F 1 (Ap A 1 (j : t))} (3.26) 



«w=2^r <fa ^ e -(( i -f)" I/, - i )-^ r (4o ■ (3 ' 27 » 

Here r(a,z) denotes the incomplete T-function. Similar expressions can be derived for 
higher odd moments \n(j, t)\ 2e+1 . 

Starting with -Fi(O) = 1, the function Fi(z) decreases monotonically, and asymptotically 
it behaves as 
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Since for t — > oo, the function Ai(j,t) ~ pot — > oo, the contribution of i*\ in this limit can 
be neglected, and we find the asymptotic result 



|nO-,OI«"?=(Mi(j,0) 1/2 = V- (^0'^)) 1/2 • ( 3 - 29 ) 

V 71 " V TT 

In the next section we will find that Ai(j,t) in a large initial region increases quite slowly 
with t, so that the simple relation (3.29) holds only for rather large time. As we will see, 
this may be one reason for the failure of previous simulations to exhibit the t 1 / 4 -behavior. 

For a Gaussian distributed continuous variable Eq. (3.29) is rigorous. The correction 
function F\(z) essentially reflects the fact that n(j,t) ranges over a discrete set of integers. 
This easily is checked by calculating |n(r)| for a process, Gaussian distributed with width 
T, but restricted to the integers. Depending on T, we find corrections to the result (3.29) 
which numerically are quite similar to the effect of Fi(z). 

For a discussion of the general motion of some segment and its interference with tube 
renewal, it is of interest to consider the correlation among n(j, i) and n(0, t). Eq. (3.2) yields 

<p 2 (n t ,n;0,j,t) = -^ <f ^ — x~ nt x~ n P({1, . . . 1, x t }, x; j, t) (3.30) 



i r 



dipt dip e~ in ^- m * ■ exp 2p (cos ip t - 1) [A^O, t) - A 3 (j, t) 



4tt 2 

+ 2p (cos ip-1) [Ai&t) - A 3 (j,t)] + 2p (cos (<p t + <p) - l)A 3 (j,t)^ . 

We here introduced 

A.(j,*)=p(^+E^a«)) , (3-3i) 

and we used Eq. (A. 7) to evaluate (^({l, • • • , 1, %t}, t). The distribution V 2 changes its 
character for t ~ j 2 . At short times t <^ j 2 , A 3 (j,t) vanishes like e~i / 4p * (cf. Eq. A. 16), and 
V 2 factorizes. 

V 2 (n t ,n;0,j,t) = A(n t ,0,t) A(n,j,t) for t/f -> . (3.32) 

As expected the two beads and j move independently. In the opposite limit of t ^> 
j 2 , A 3 (j,t) diverges proportional to t (cf. Eq. A. 12), whereas the combinations Ai(j',t) — 
A 3 (j,t) with j' = or j, saturate at some j-dependent value. V 2 factorizes into a Gaussian 
distribution of n + = (n t + n)/2 and a distribution of n_ = n t — n: 

V 2 (n t ,n;0,j,t) = Unp A 3 (j,t))~ 1/2 exp (- - 2± — - ) Vt\n_, j, t) for t > f . 
v ' \ A Po A 3 (j,t)J 

(3.33) 

takes the form characteristic for a single bead process (cf. Eq. 3.22). 

V[-\n_,j,t) = e- A ^)/| n _|(A(j,t)) (3.34) 
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j for t ^> j z 



(3.35) 



The result (3.33)-(3.35) expresses the fact that most of the defects created by chain end 
and not destroyed again at this end by time t, for t ^> j 2 have diffused over bead j without 
returning, so that the stochastic variable n(j, t) then closely follows n(0, t). For t ^> iV 2 it in 
particular implies that the total chain motion is dominated by the single stochastic process 
n(0,t). This observation underlies all theories of the large time limit f[|,l4] . 



IV. MOTION INSIDE THE TUBE 



Within our model, we now can treat the thermal displacement of some bead j exactly as 
long as t is so small that the tube renewal has not reached that bead. Eqs. (2.8) and (3.26) 
yield 



g 1 (j,N,t) = ((r J (t)-r ] (0)y) 



e.^n*(j,t) (l-F 1 (2n 2 (j,t))) 



(4.1) 



n 2 ( J} t) =2p A l (j,t) 



(4.2) 



In view of the symmetry of the chain we henceforth restrict the analysis to < j < y. We 
now evaluate several limits of g\ (j, N, t) analytically. 

A) Introducing the variables 



pt 
N 2 ' 



J=J/N 



(4.3) 



we take the limit of long chains iV — > oo, and we consider intermediate times < 
i < oo, i.e., times of order T 2 . We thus focus on the regime T C f < T 3 , where 
the crossover from to g\ ~ t 1 / 2 occurs. (See Eq. (1.2).) It will be shown in 

Sec. V.C, that for j > it needs a time i = Tr(j') — > oo, before the tube is destroyed 
up to bead j = jN — > oo. Thus in the present limit tube renewal is irrelevant and our 
results are exact. 



It now is useful to rewrite Ai(j,t) as (cf. Appendix A.l) 

A 1 {j,t) = NA 1 Q,t) 



(4.4) 




4p sin2 



(4.5) 



cos 2 (irkj + ^) 

^ (i) 
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Taking N — > oo, we find the limiting form of Ai(j,t): 



l t r 1 ds 



+00 



2U Jo sfs v ^ 



„2 



e ^ + e 



(4.7) 



Eq. (4.7) results from a Poisson transformation and is easily expressed in terms of 
error functions. 

Clearly Ai(J,t) diverges for iV — > 00, j,i fixed. Therefore the correction function 
Fi(. . .) in (4.1) vanishes, and we find 

g 1 (J,t) = lim N- 1 ' 2 g 1 {j 1 N 1 t) 

N_ — » 00 
t,j fixed 

= 24^A(j,t)) 1/2 (4.8) 

For chains of length N finite, but so large that they can be described by these limiting 
expressions, we obviously can absorb all microscopic model parameters into a rescaling 
of N and t 



2 

s 

t^t = £jp 2 pt (4.9) 



N -> N = t s p N 

(A J. 



Thus the present results are independent of the microstructure and can be written in 
universal form as 

9l (j,N,t)=^%(l,^j . (4.10) 
They also can be derived from a continuous model and agree with the result of ref. 

1 /2 

[|], appendix B, except that in that work gi was identified with (n 2 (j,t)^j , in our 



language. This introduces an additional factor yvr/2. 

Considering now extreme values of i we from Eqs. (4.6), (4.7) find the power laws 



t ,t > 1 



Note that the limits i — > 0, j — > do not commute. Taking j 1 — > first, we from Eq. 
(4.7) find i.e., twice the result of the first line in Eq. (4.11). Below we will see 

that this induces additional structure in the t 1 / 4 - regime. Eqs. (4.8), (4.11) yield 
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- gi Q,t) = n s M.\ i") V4 rt<<1 ' (4.i2) 



7T 



t 1 / 2 ,t»l 



equivalent to the well known asymptotic results 



1/4 

9l (j,N,t)~-^ ■{ y> ' ■ (4.13) 



2_ (1)- , *<1 

V^F | #-1/2 £1/2 ^ f>1 



For any j > 0, the limit t — ^ is independent of N, since for t <^ T 2 — 0(N 2 ) bead 
j = jN does not feel the presence of the chain ends. Its motion is driven by the 
particles which initially sit in its neighbourhood on the chain. In contrast, for t ^> T 2 , 
the motion is dominated by particles diffusing over the total length of the chain, which 
yields the factor iV -1 / 2 , again independent of the position of bead j inside the chain. 

The crossover among the limits of small or large time depends on j. Naively we might 
expect it to be governed by the time needed for particle diffusion from the nearest 
reservoir to bead j: 

f 2 Q)~3\ ®<J<\ . (4.14) 
For < j <C \ the situation, however, is even richer. The crossover 

£l/4 ^ (1/2 a l ways 

occurs at 

t wf 2 Q wO(l) . 

It is due to particles diffusing over the whole length of the chain with nonnegligible 
probability. For such times several terms nonnegligibly contribute to the sums in Eqs. 
(4.6), (4.7). For j <C |,T 2 (j) ^ t < f 2 Q), however, the v — contribution of the 
second term in Eq. (4.7) is of the same order of magnitude as the v — contribution 
of the first term, both being much larger than all v ^ terms. Thus the t 1 / 4 regime 
splits according to 

*<3.i> = U*l T [-) f 2 (])«f«f 2 (l) 

The crossover at T 2 (j) is due to the particles diffusing in from the nearest chain end, 
while the other chain end effectively is still infinitely far away. 

The sums in Eqs. (4.6), (4.7) are rapidly converging in overlapping domains of i. gi(j, t) 
therefore is easily evaluated numerically. Fig. 3 shows typical results, normalized to 
the asymptotic t — > behavior. The crossover from g\ ~ t 1 / 4 to g~\ ~ t 1 / 2 occurs for 
0.1 < i < 1. For j = 0.01 the splitting of the t 1 / 4 -regime is clearly seen. 
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B) We next consider the limit of long chains iV — > oo, intermediate beads < j < 1/2 and 
shorter times < t < oo. In this limit for beads j = jN ^> 1 both chain ends effectively 
are infinitely far away and the result becomes independent of j. We thus describe the 
behavior of some bead well inside the chain for times t <C T 2 ~ 0(N 2 ) — > oo. The 
analysis, based on Eq. (3.12), reveals that A\(j,t) reduces to (cf. Eq. (A. 5)) 

Ai(j,t) = -pt dx cos 2 x(l -4p sin 2 a;)*" 1 , (4.16) 



with the asymptotic limit 

1/2 



Al (j,t)~(^) , t»l . (4.17) 



In this limit we recover the results found above for t < 1. 

The dominant corrections to that behavior arise from the contribution Fi(ApoAi(j,t)) 
in Eq. (4.1). With Eqs. (4.17), (3.28) we find 



9i(j,N,t) 



24 (- 

7T 



Po\ 1/2 fpt\ 1/4 



71 ) 



1 



16p Vpw 



(4.18) 



With p ~ 1/4 the leading order correction in (4.18) falls below the level of 1% only for 
pt > 5 • 10 3 . Fig. 4 (full line) shows our full result, normalized to the asymptotic t 1 / 4 
behavior. We choose values t s — 2,p = |, p — j, but we should note that for given 
i s the result essentially depends on the scaled variable p\pt ~ t, only. The residual 
Po-dependence is negligible for p\pt > I. This reflects the fact that A 1 (Eq. 4.16) 
rapidly attains its asymptotic behavior (4.17) the corrections dropping to the level of 

1 % for pt ^ 10. By virtue of Eq. (4.2) this implies that (n 2 (j : t)j rapidly tends to 
its asymptotic t 1 / 4 -behavior. The variation seen in Fig. 4 solely is due to the correction 
function Fi(z), that reflects the discreteness of n(j,t), cf. discussion after Eq. (3.26). 

This long initial transient, extending almost up to pt = O(10 4 ), at a first glance is 
quite surprising. A second thought reveals that it is quite natural. With the estimate 
\n(j,t)\ ~ (pt) 1 / 4 we find that in the region of interest typical values of n(j,t) are of 
the order 10 or smaller. It is clear that with such small values of n(j, t) the discreteness 
of the process may induce large corrections. 

We now also might discuss the limit iV — > oo and j, t fixed, implying that the bead considered 
has a finite distance from a chain end. Here, however, tube renewal rapidly comes into play, 
and we therefore postpone the analysis to the next section. We close this discussion by 
illustrating the motion of the central bead j = N/2 for times t C T3 in a chain of finite 
length. The results follow from numerical evaluation of Eqs. (4.1), (4.2), with Eqs. (3.27), 
(4.4), (4.5) inserted. Fig. 4 (broken lines) shows the typical behavior of g\ (^,N,tj, again 

normalized to the intermediate asymptotic t 1 / 4 -behavior. For iV = 10 3 this behavior is 
attained over almost three decades, until for pt > 10 5 the curve bends over towards the 
t 1 / 2 -regime. With decreasing iV the length of the t^-plateau decreases rapidly, and is just 
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about visible for N = 10 2 . For shorter chains it is washed out. In a doubly logarithmic plot 
the result for N = 40, for instance, for short times would look like some effective power law 
close to gi ~ t 1 ^ 3 . 

As shown in Fig. 3 the crossover towards the t 1//2 -law starts for pt « iV 2 /10 and takes 
more than one decade in t. The t 1//2 -behavior reasonably well is approached for pt ^ 2N 2 . 
In the next section we will find an estimate for the reptation time: pT% « ^ N 3 . Thus the 
condition for the occurence of a true t 1//2 -regime: T2 <C t <C T3 takes the more quantitative 
form 2N 2 < pt <C j= N 3 . Clearly for typical chain lengths used in simulations (N < 10 3 ) 
this condition is not fulfilled in a sufficiently large range of pt to allow for an intermediate 
t 1 / 2 -region. 

Above we stressed that the initial slow transient in gi(j,N,t) is due to the correction 

Fi(. . .). (n 2 (j,t)^j accurately shows t 1//4 -behavior starting as early as p^pt = 0(1). This 
quantity therefore shows standard reptation behavior much clearer than gi(j, N, t), and it is 
of interest to ask whether it can be related to some simple observable. Indeed it is not hard 
to show that for the lattice model described in Sect. 2, it is related to the cubic invariant 

\a=l 

= i s (^U7t)) 1/ \ t<T R (j) . 

This result holds in three dimensions. For other dimensions or other microstructures 
(n 2 (j,t)^j always can be measured as an appropriate combination of the fourth and the 
second moment of (rj(t) — Tj(0)). 



V. TUBE RENEWAL 

Tube renewal is just another word for the motion of the end beads that for longer times 
also affects the inner beads and the chain as a whole. As shown in Sect. II. A for chain end 
and times t < T3 it leads us to consider n max (t), defined as the maximal value within 
time interval [0,t] of the stochastic variable n(0,t), that gives the change in the occupation 
number of reservoir (cf. Eq. (2.10)). In particular the relation (2.12) 



gi {Q,N,t)= (r (t) - r (0)) 2 



2£ s n max (t) 

holds. 

As is obvious from Eq. (3.13), n(0, s) executes a correlated random walk with n(0, 0) = 0, 
the correlation arising from the contribution ~ ^2(^2 — Si): a particle, that jumped out of 
the reservoir at time s±, with some probability depending on S2 — s\ > will fall back again 
at time S2- For correlated processes no manegeable exact expression for n max (t) is known. 
Only for uncorrelated walks a result known as Spitzer's identity |HJ yields a simple answer. 



(This identity is further discussed in appendix B.) For a reflection symmetric uncorrelated 
walk n s starting at no = 0, it yields the relation 
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rimaxit) = £ if- 
s=l ZS 

We will base our analysis on this result. To get some feeling for the structure of the problem 
and the approximation involved we, however, first derive a formally exact expression for 
n m ax{t)- We then approximately treat the motion of a chain end, and with a further ap- 
proximation we calculate the motion of some bead j for times t > Tr(j'). Finally we discuss 
the reptation time. 



A. Formal analysis of n max (t) 

From Eqs. (3.2), (3.17), (3.13) the distribution of {n s } = {n(0,s)} is found as 

P 2 (K};0 = (27T0-* / n ^x: n ° expQ 2 ({s.};f) , (5.2) 

J s =l Xs 

where all integrals can be taken to range over the unit circle. It is useful to transform from 
variables x s conjugate to n s to variables y s conjugate to the increments n s — 

x s = 1 < s < t; y t+1 = 1 (5.3) 

Vs 

Successive substitution yields 

V 2 ({n s };t) = (27T1)-' <fti — (—) S expQ 2 ({y s };t) , (5.4) 

j s =i y g Ws+ij 

where 

Q2{{y s }] t) = Po pJ2 (y s + --2) 

s=l V Vs J 

+PoP 2 E ( — - 1 ) (y-2+1 - !) M*2 - si) (5.5) 

0<S!<S 2 <t \^ s l / 

To derive an expression for n max {t) we define the operator 

t 

B(n,t) = [] e(n-n 8 ) , (5.6) 

s=0 

with the discrete function as defined in (2.14). B(n,t) obeys the relation 

B(n,t)-B(n-l,t)=6 nnmax{t) . (5.7) 

We thus have 



,(*) = £ n (5(71,0-5(71-1,^) . (5.8) 



n=l 
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We next introduce 



C(t) = n max (t) - n max (t - 1) , (5.9) 

so that 



n max (t) = J2 C(s) . (5.10) 

s=l 

Combining Eqs. (5.8), (5.9) we easily find 

oo 

C(t) = -J2 (WT)-B(n,t-1)) 

n=l 

oo n oo 

= E E E Min.Y,t) • (s.ii) 

n=0 ni...nt-i=— oo nt=n+l 

In Eq. (5.4) for V2({n s };t) we now order the integration contour B as 

B: 1 > |yi| > |s/2 1 > ••• > \vt\ (5.12) 
We then in Eq. (5.11) can carry through all summations to find 

f n _^ l_ 



C(*) = (27T1)-* £ n - ^ ~~T ex P g 2 ({y s };t) (5.13) 



Recall the definition = 1. 

Eq. (5.13) is an exact but essentially intractable expression. To establish the connection 
to Eq. (5.1) we note the identity 

(27ri)-* / n — ^ J- £ 5 s ,F{y s } - -L / rf y F | y ' • ' ' ' y} = , (5.14) 

JB s =i y s - y s +i ^ 2tt? /6<| w |<i (y - l) 2 

where <S S / is a cyclic shift operator: 

S s >F{y 1 ,...,y t } = F{y 1+S ,,y 2+S >,. .., y t ,Vi,y2, ■ ■ ■ ,Vs>} • (5.15) 

This identity holds for any function F{. . .} analytic in all variables y s in the domain 

a > \y s \ > b ; a > 1 > \y t \ > b 

We now note that 

Q 2 ({y,...,y};t) = Q 1 (y;0,t) , 

which yields 

1 r expQ 2 ({y,...,y;t) 1 



, . t dy y T i\ 2 =il"(°»01 (5-16) 

2?n /|„|<i (y - l) 2 2 
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(cf. Eq. 3.25). We then subtract from C(t) (Eq. 5.13) the identity (5.14) multiplied by 1/t, 
to find 



e{t) = MM 



+{2m)-* I f[ — — r fl " 7 £ ^ ] expQ 2 ({2, s };t) (5.17) 




We note that the contribution \n(0,t)\/2t takes into account the strongest singularity of the 
integrand in Eq. (5.13), occuring if all variables y s coincide. 

For an uncorrelated random walk the exponential factorizes: exp Q2{{ys}, t) — > 
ris=i f(Vs), and the second contribution in Eq. (5.17) vanishes. In view of Eq. (5.10) 
the result therefore yields a formal generalization of Spitzer's identity (5.1) to correlated 
walks. Transformed back to the variables n s , the second term in Eq. (5.17) weights each 
walk contributing to C(t) (Eq. 5.11) by the difference of its proper weight and the average 
weight of all of its cyclic permutations. Depending on the walk considered, the result can 
be of either sign, and we may hope that cancellations keep the contribution of the second 
term in Eq. (5.17) small as compared to the first term. We, however, have found no proper 
estimate. 

More can be said only in the limit iV — > oo,t = pt/N 2 fixed (case A of Sect. 4). For 
s = ps/N 2 ^> 1 the coupling term ^(s) in Q2 (Eq. (3.13)) vanishes exponentially, reflecting 
the fact that over time intervals of order T 2 ~ iV 2 the particle distribution within the chain 
equilibrates. A process with coarse grained time intervals r ^> T 2 proceeds by steps of 
average size (An) 2 ~ r ^> 1, with correlations of order 1 of a time step with its immediate 
precursor. Neglecting these small correlations we find that Eq. (5.1) gives the leading con- 
tribution to n max (t) for i — ► 00. In this long time regime chain motions can be modelled as 
an uncorrelated random walk of a single degree of freedom 11,14]. 



B. Motion of the end bead in the 'mean hopping rate' approximation 

We here evaluate <?i(0, N, t) as 

gi(0,N,t)=£, J2 -RMI , (5.18) 
s =i s 

(cf. Eqs. (2.12), (5.1)), with 

|n(0,s)| = ^= (fioA^O, s)) 1/2 [1 - F 1 (4p A 1 (0, s))} , (5.19) 

(cf. Eq. (3.26)). 

Since in Eq. (5.18) the contribution C(s) of a time step s is calculated as the contribution 
of an uncorrelated walk giving the same (time dependent) distribution of n(0, s) as the true 
correlated walk, we call this a 'mean hopping rate' approximation. For the interpretation 
of our results it is useful to determine the effective hopping rate p e ff(s) of this uncorrelated 
walk n e ff(s). Writing 
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n 2 eff (s) = 2 Peff s 



and identifying n 2 e ^{s) and n 2 (0, s) we find 

p eff (s) = ^A 1 (0,s) . (5.20) 

Eq. (3.12) yields A±(0, 1) = p, so that the expected result 

p e ff(s = 1) = pop (5.21i) 
follows. For large s, however, the first term in Eq. (3.12) dominates and we find 

p eff (s -> oo) -> — (5.21n) 

Thus the effective hopping rate strongly decreases with increasing time, reflecting the fact 
that particles emitted from the reservoir may have been absorbed again. 

We now first consider the motion of the end bead (5.18), (5.19) in several limits, corre- 
sponding to those of Section IV for some internal bead. 

A) iV -> 00, t fixe d. 

Since |n(0, s)\ increases with s, the sum (5.18) is dominated by large values of s. In the 
limit considered, we therefore can neglect the correction Fi(. . .), and we can replace 
the summation by an integral. This yields 

gi(0,i)= lim N-^ 2 9l (0,N,t) 

N -» 00 
t fixed 

= U.\I^ t-^^S) , (5.22) 

where 



1 Jo s 



, 2J± , s < 1 
Ax(0,s)~^ V?r (5.23) 
s , s»l. 



(Cf. Eq. (4.11) and the associated discussion.) 
For [> 1 we find 



gi(0,I) = U.(—)' , (5-24) 



7T 



a result which in the light of the discussion of the previous section is expected to be 
exact. Comparing to the motion of an inner segment j > 0, Eq. (4.12), we find 



gi(M) *^?o 
<7i(J,*) 



2 . (5.25) 
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Recall that due to tube renewal the effective random walk of an end bead in the limit 
considered has twice the length of the walk of an interior bead, cf. Eqs. (2.8) and 
(2.12) with n max {t) ~ \n(j,t)\. This is the origin of the result (5.25). 

For t < 1 Eqs. (5.22), (5.23) yield 

^ 1 (0,t)=8v / 24f^V /2 f-V » ( 5 - 26 ) 



7T / V 7T 



or 



gi (0,*) 



'=4V2 . (5.27) 



This result is not trivial. Naively we would estimate the ratio (5.27) as 2 a/2, the 
factor 2 reflecting tube renewal and the factor \/2 arising from the enhanced mobility 
of the end bead. (Cf. the Poisson representation (4.7) for t — > 0,j — 0, as compared 
to j > 0.) The additional factor of 2 found in Eq. (5.27) is due to the integration over 
s (Eq. 5.22), reflecting the fact that n max depends on the history of the walk. In our 
approximation this enhancement of the ratio (5.27) can be seen as a direct consequence 
of the slow decrease of the mean hopping rate discussed above. 

Fig. 5 shows the crossover behavior of the ratio g~i(0, t)/g~\ (j},t). The crossover from 
the asymptotics (5.27) to (5.25) sets in at t 0.1, i.e., t T 2 /10. It is quite slow, 
taking about three decades of t. We should note that at time t ~ T 3 , formally infinite in 
the present limit, there must be another crossover towards gi (0,N,t)/ gi (%,N,t) = 1, 
since eventually the chain is displaced as a whole. 

B) iV -> oo, t fixed 

This limit shows the approach towards the asymptotic law (5.26), now adequately 
written as 

01(0,00,*)^ 8V2i s ^) 1/2 (^) 1/4 . (5.28) 

In the general expressions no special simplifications occur, except that A 1 (0,t) (Eq. 
(3.12)) takes the form 

/ \t— 1 

Ai(0,t) = 8pt / dx cos 2 (7ra;) [1 - Ap sin 2 (7ra;)J 
fl/2 \t 1 

+ J dx [l-Ap sm 2 (7ix)j -- . (5.29) 

The result for gi(0, oo, t), normalized to the asymptotics (5.28), is shown as full line in 
Fig. 6. As compared to the behavior of an interior segment (Fig. 4) we find a drastic 
enhancement of the initial effects. (Note the change in the scales.) The t 1/,4 -regime is 
reached only for pt > 10 8 . This again is due to the discrete nature of the process, which 
here comes in at two places: Firstly, \n(0, s)\ involves the correction F 1 (4p A 1 (0, t)) 
(cf. Eq. 5.19), and secondly, n max (t) involves the sum over all time steps. The long 
initial transient is dominated by the latter feature and thus quite insensitive to the 
value of po, that only influences Fi(. . .). 
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Evaluating gi(0,N,t) for finite N and normalizing by the intermediate asymptotics we 
find results as shown by the broken curves in Fig. 6 . For the end bead it needs extremely 
long chains to see indications of t ^-behavior. For chain lengths N < 100, typically used 
in computer experiments, the initial transient directly bends over towards free diffusion. 
Clearly the observation of an intermediate t 1//2 -regime is beyond all reach. 

C. Motion of an arbitrary bead 

As shown by Eq. (2.13), the motion of some bead < j < N for times t < T3 is a 
superposition of motion within the tube and tube renewal. This leads to a rich structure, 
which we first discuss on a qualitative level. A schematic plot, summarizing our discussion 
for some bead 1 <C j <C N/2, is shown in Fig. 7. 

Consider some bead 1 <C j ; <S y in an extremely long chain. Initially it does not feel the 
chain ends but shows the slow transient towards a t 1//4 -behavior, as discussed in Sect. 4: 

9l (j,N,t)^2£ s M fPLY f \ T «t«T 2 (j), pT 2 (j)^f . (5.30) 

V 71 \ 71 / 

For t ~ T 2 (j), particles coming from the nearest chain end start to diffuse over bead j, which 
raises the amplitude of the t^-law by a factor of \[2, cf. Eq. (4.15), thus 

9l (j,N,t) *2V2£ S M f^V /4 , T 2 (j) «t «T R (j) (5.31) 



7T \ 7T / 

This regime will end at a time T R (j) defined as the time needed for the tube renewal to pass 
over bead j: 



isn max (T R (j))=j (5.32) 
On the basis of Eq. (5.28) T R (j) is estimated as 

For t > Tr(j'), tube renewal dominates. For extremely long chains it yields another t 1//4 - 
regime. 

9iU,N,t)^8V2£ 8 M f^Y /4 , T R (j) «t «T 2 (iV) (5.34) 

V 7T \ 7T / 

Gradually the t 1 / 2 -behavior sets in beyond the Rouse time T 2 (N): 

9l (j,N,t)^A£ s (~) 1/2 , T 2 (iV)«t«T 3 , (5.35) 
and finally we end up with free diffusion 

9l (j, N, t) w const ^, T 3 < t (5.36) 
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T 3 may be identified with Tjj(y), calculated from Eq. (5.32). For j = O(y) we, however, 
can not use Eq. (5.33). In evaluating Eq. (5.32) we rather have to resort to Eq. (5.24) to 
find 

7r i 2 N N 

The above analysis breaks down for Tr(j') > T 2 (A r ). With T 2 (N) ~ N 2 ,T R (j) ~ j 4 
we find it to be valid only for 1 <C j <C Oi^N 1 / 2 ). For larger values of j the t 1//4 -regime 
(5.34), driven by the motion of the end bead, vanishes. Then tube renewal under favourable 
conditions may introduce some structure in the t 1//2 -regime. We, however, should note that 
this whole discussion is more a matter of principle. It will need extremely long chains 
(N > 10 4 ) to observe these different regions. In practice we expect that only the initial t 1 / 4 
regime (5.30) properly can be observed, the other regimes being covered by broad crossover 
regions. 

For a more quantitative analysis of g±(j, N, t) as given by Eq. (2.13) we have to construct 
the simultaneous distribution of n max {t) and n(j,t). 



'P m ax,j(n m ,n;t) — 5 nmax ( t ),n m ^n(j,t),n (5.38) 

Again the formal expression resulting from Eqs. (3.11) - (3.17) is intractable. To construct 
an approximate form we introduce n t = n(0, t) as intermediate variable and write 

+00 

'Pmax,j{n m ,n;t) = ^ ^n max ,n m ^n(0,t),n t ^n(j,t),n 
nt=—oo 
+00 

= 'P m ax,o,j(n m ,n t ,n;t) . (5.39) 

nt=— 00 

We then factorize according to 

j(n m ,n t ,n;t) ^V ma x,o( n m\nt;t)V 2 (nt,n;0,j,t) , (5.40) 

where Viijit-, n; 0, j, t) is the simultaneous distribution of n(0,t),n(j,t) discussed in Sect. 3, 
and Pmax,o( n m\nt',t) is the conditional probability to find n max {t) = n m , given n(0,t) = n t . 
This approximation serves to isolate the hard part of the problem in the analysis of the 
single stochastic variable n(0,s). It is justified by the observation that for beads not too 
close to a chain end (j > 20, say), the relation T R (j) ^> T 2 (j) holds, so that bead j feels 
the tube renewal only for times where n(j,t) is bound to n(0,t) fairly rigidly. (Recall the 
discussion after Eq. (3.35).) 

To find an approximate expression for the conditional probability V maXt o(n m \n t ; t) we 
resort to the simple random walk, which is the only model where this probability can be 
calculated in explicit closed form. We thus ignore all the correlations of the stochastic 
process n(0, s), but we keep the important inequalities < n max (t), n(0, t) < n max (t), and 
we can adjust the hopping probability of the random walk so as to reach the desired values of 
n max {t). As a final purely technical point we take the variables n m , n t , n to be of continuous 
range, which is a minor approximation since for beads not too close to a chain end we are 
interested only in fairly large times {pt > pT R (j) > 10 3 ). 



26 



With these approximations the explicit construction of V max j(n m , n; t) is carried through 
in appendix C. The final result reads 



e( Zrn ) 



(t) (n 2 C/,t)) 



1 

V2 - 



2U/2 



exp 









I -a 2 
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+a ( 2a\/- z m — z\ exp 



where 



-- I 2aW- -2m - -2 

7T 



erfc 



7T ^ 



'l-2a 2 ) + az 



^(1 " « 2 



(5.41) 



7l T , 



:(*)' 



and 



a = a(j,t) = 



(n 2 (j\0) 1/2 ' 



(^(o^AO',*)) 1 / 2 



(5.42) 



^(j, t) is defined in Eq. (3.31). The parameter a, < a < 1, measures the strength of the 
correlation among n(0,t) and n(j,t). For a — > 0, which is reached for t <C T 2 (j), the two 
processes are independent. They are tightly bound together for a — > 1, corresponding to 
f>T 2 {j). 

To give an impression of the shape of the distribution we in Fig. 8 have plotted 
V(z m ,z,t) = n rnax (t)(n 2 (j,t)) l/2 V maxJ (n m ,n;t) for fixed z m = 1,2,3 as function of z, 
using values a = 0.9,0.99. Specifically for t = T R (j), i.e., for times where the tube re- 
newal just reaches segment j, such values a = a(j,T R (j)) correspond to j = 50 or 120, 
respectively. As is shown, with increasing a the curves become strongly asymmetric. This 
reflects the distribution of n(0,t), which is cut off sharply at n(0, t) = n max (t). We also 
note that tube renewal will have a substantial effect even for t < T 3 (j). Both conditions 
is n max (t) > j, i s n max (t) > j + i s n(j,t) then are fulfilled for z m > l,z < 0, a region that 
carries a considerable fraction of the weight of the distribution V max ,j- 

We now evaluate gx(j,N,t), decomposing Eq. (2.13) as 



9l (j,N,t) = 9i (j,N,t) + g r (j,N,t) , 



where 



gi (j,N,t)=i s \n(j,t)\ 



(5.43) 
(5.44) 



is the contribution of the motion in the fixed initial tube, as discussed in Sect. 4. The 
contribution of the tube renewal can be written as 



g r (j,N,t) = 2£ s e(-n(j,t) - 1)0 (n rnax (t) - j-^j (n rnax {t) - j-^j 



+ 2i s e(n(j, t))e (n max (t) - n(j, t) - j-^j (n max (t) - n(j, t) - j-^j (5.45) 
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With our appoximations it yields 



9r(j,N,t) = 21 s | (t)\ J dz (t)\ ,0,a(0,j',t) 



+G 



(*)l > ■ ,L 



l^max (^) | 



J 



(5.46) 



G(z, 6, a) 



- (l-2ab + b 2 ) - 

7T V 7 1 



1/2 1 — afe 



• erfc 



- 2ab 
b — a 



cxp 



- — fl -2ab + b 2 

71 V 



-1 



2" 



(1 - 2abY x exp 

7T 



-4a' 



,2 



7T 



erfc 



I -2a 2 



/tt(1 - a 2 ) 



(5.47) 



Typical numerical results are shown in Fig. 9. Together with the full results for iV = 10 3 , j = 
20 or 100, we included . .) and g r (. . .), and we show <?i(0, iV, t), for comparison. In . .) 
the first splitting of the t 1//4 -regime can be seen both for j = 20 and j = 100, but g r (. . .) 
only for j = 20 shows weak indications of the t 1 ^-behavior from the end bead motion. 
For j = 100 tube renewal is effective only for t > T 2 (y) . In both cases we see that the 

total motion g\(j,N,t) beyond the first t 1//4 -regime is dominated by the crossover. We also 
note that tube renewal effects set in for t > T R (j)/ 10, a consequence of the fairly broad 



distribution V, 



max,] 



D. An estimate of the reptation time 

The reptation time T3 measures the time needed to destroy the original tube. This process 
proceeds by the motion of both chain ends, and a simple definition consists in identifying 
T 3 with T R (N/2), Eq. (5.31): 

N 

4 n max (T 3 ) = — . (5.48) 

Other, more complicated definitions might be given, that take into account the first-passage- 
time nature of the problem. The present definition closely is related to the evaluation of the 
melt viscosity in the reptation model [rjj,n^|T3]. With Eq. (2.12) it implies 



9l (0,N,T 3 ) = N = R 2 e , (5.49) 

where R 2 is the mean squared end-to-end distance of the chain. We should note, however, 
that strictly speaking for t = T3 Eq. (2.12) is invalid due to the interference of chain renewal 
from both ends, so that the physical interpretation of T 3 , as implied by Eq. (5.49), is only 
approximate. 

For long chains we may invoke Eq. (5.37) to find the well known asymptotic iV 3 -law: 
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(5.50) 



To estimate the corrections we first may evaluate Eq. (5.48) with the form (5.22) of 
(- s n m ax{t) = iV 1//2 (7i(0, t)/2, valid for long chains. This is equivalent to Doi's model fl5| 
and yields the thin line in Fig. 10. As is seen this curve bends downwards, yielding a strong 
increase in the effective exponent T 3 ~ N x with decreasing N. In fact, naively extending 
this curve towards N — > 0, we approach x = 4, reflecting the t 1//4 -law. Of course this exten- 
sion is meaningless since T^(N) by construction is far beyond the range where the t 1//4 law 
is valid. It is not surprising that the result is strongly modified if we take the discreteness 
of the system into account. Our result using the full expression 5.1 for n max {t) is shown 
as the fat line in Fig. 10. For N > 200 it merges with the continuum approximation, but 
for the full range 20 < iV < 1300 it is nicely approximated by pT 3 (N) « 0.059iV 3 - 15 . We 
thus see that the slowing down of the increase of n max (t) with increasing t implied by the 
correlations increases the effective exponent beyond x = 3, but not as strongly as expected 
from a continuum approximation. 

Evaluating the Doi-Edwards expression for the melt viscosity rj we find results for rj(N) 
most similar in all respects to those shown for T^(N) in Fig. 10. For 20 ^ iV < 10 3 the 
full theory yields an effective power law rj(N) ~ N 3 ' 2 , with the continuum approximation 



bending downwards as in Fig. 10. Note that these results are quite similar to those of |]18 
Fig. 3, where Doi's 'Rouse chain in a tube' model was simulated. 



VI. CONCLUSIONS 

In this work we evaluated a discrete version of the reptation model, concentrating on the 
motion of the individual beads. Depending on the bead position along the chain we found a 
very rich structure, characterized by fairly broad successive crossover regimes. Initial effects 
due to the discreteness of the model die out very slowly. For chains of length N ^ 100 they 
are visible for the central bead up to times of the order of the Rouse time, and for the end 
beads they are even more pronounced. 

Our model assumes the smallest possible tube diameter in the sense that chain segments 
between two successive entanglement points show no internal structure. Even in the nonuni- 
versal initial range the results therefore should be relevant for the quantitative interpretation 
of simulations of lattice chains in surroundings of high obstacle density. Indeed, as shown 
H and in our subsequent paper |ll|, they perform very well in explaining data of the 



m 



Evans-Edwards model. A glance to Fig. 8 of ref. || or Fig. 6 of ref. |T{| reveals that also 
more complicated dense lattice or continuum systems for the motion of the central bead of 
the chain yield results closely resembling our findings. We therefore believe that even in 
the microstructure dependent range our results quite generally show the typical effects of 
reptational motion. 

There remains the question which features could be used to distinguish reptation from 
motion dominated by disorder or affected by relaxation of the environment. As we have 
seen, the unambiguous identification of the power laws will be most difficult. For values 
N = O(100), as typically reached in simulations, only the motion of the central bead shows 
a small region described by the t 1//4 -law, and to unambigously identify this behavior as a 
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limiting law its stability with increasing chain length would have to be checked. Beads not 
close to the center of the chain feel successive crossovers which essentially wash out all the 
t 1 / 4 -structure. However, as pointed out at the end of section 4 and illustrated in fTT| , |T2f . 
appropriate fourth moments of Tj(t) — i"j(0) show the t 1//4 behavior much more clearly, even 
for very short chains. Such moments therefore seem to be the appropriate objects in a search 
for the initial t 1//4 behavior. At least, a comparison of the second and fourth moments will 
allow us to identify the range and magnitude of the nonuniversal initial effects. 

According to reptation theory the t 1//4 -regime for T 2 C t C T 3 is followed by t 1 ! 2 - 
behavior. Our results show that for typical chain lengths this interval is too small to 
identify such behavior, in particular if we take into account the fairly strong fluctua- 
tions in the tube renewal (n max (t)). In practice we for t > T 2 expect to see a smooth 
crossover towards free diffusion. To find indications of the i 1//2 -regime it will be preferable 
to measure the motion of the central bead relative to the center of mass: g2(^ ) N 1 t) = 

((r N / 2 (t) — R cm (t) — rjv/ 2 (0) + R cm (0)) 2 y This quantity for t <T 3 shows the same behav- 
ior as gi(^, N,t), but an eventual increase of the effective exponent from x ~ 1/4 towards 
x = 1/2 cannot be interpreted as crossover towards x — 1, since g 2 (- ■ ■) for t > T 3 crosses 
over to a constant g 2 {^ 1 N,t) « R 2 g . 

As a further quite pronounced signature of reptation we note the difference in the be- 
havior of the end bead as compared to the central bead of the chain. For t -C T 2 the 
ratio gi(0, N, t)/(?i(y, N, t) in the Rouse model tends to 2, and we do not expect this to be 
strongly modified by disorder, which in leading approximation just renormalizes the over all 
time scale 0,0- In contrast, we in our mean hopping rate approximation find a much larger 
asymptotic value 4 a/2 of this ratio. This not only illustrates the effectiveness of the tube in 
constraining the motion, but also it is a consequence of the slow decay of the correlations 
of n max (t) in time implied by the basic defect motion. We believe that this enhancement of 
<7i(0, N, t)/gi(y, N, t) should be observable fairly easily. 
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APPENDIX A: SOME PROPERTIES OF THE COEFFICIENTS A 



Eq. (3.12) defines 

Evaluating the first contribution in the sum we find the equivalent form 

, , ■ n Pt N 1 1 / 1 \ . j 2 

A 1 (j,t) = —+ + 1 ? + — 

u ' N 3 2 6N V N J N 



1 w 



Defining Al = Ai(j,t)/N; j = j/N,i = pt/N 2 we can take the continuous chain limit 
N — ► oo, i,j fixed, to immediately arrive at the results 



7T JO 



e *» + e 



(A4) 



(For comparable calculations in more detail, see the appendices of J2IJ.) Besides yielding 
the asymptotic limits £ — > 0, oo, these expressions allow for an efficient numerical evaluation. 
Up to errors less than 10~ 15 we may use the first form truncated at k = 3 for i > 0.2, and 
the second form truncated at \v\ — 2 for t < 0.2. 

In the limit N — > oo, t, j fixed, we go back to the form A.l written as 

N-l 



1 + COS ( 7T k ( 2j + 



1 

N 



2^ Jo 



*/* , 1 - (1 

ax 



4psin x) 



sin 2 x 



1 ^ 1 



+ 4AT . 



<4 



fe=i sin ^2]v 



COS ( 7T fc ( 2j + — 



The last contribution for N — ► oo oszillates to zero. The first term by partial integration 
can be written as 



4 / >7r / 2 / 
Ai (j; t) = — pt dx cos 2 x ( 1 — 4p si 
n Jo v 



sin 2 x 



t-i 



(A5) 
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2. A 2 (s) 



Eq. (3.14) defines 



N-l 



Ms) = tfE 



, , at sin 2 



~N 



It is straightforward to derive the relation 



E Ms) = I 



1 - 



1 



JV-1 



NJ 2p 2 N ^ 



E (i - «*) 



cos 



2 I nk 
2N 



sin' 



leading to 



pt-p 2 E A 2 (s) = ^!(0,t) 

0<si<s 2 <t 



(A6) 



(AT) 



For s > iV 2 , A 2 (s) vanishes as exp (-tt 2 ^)- In the limit N — > oo, s fixed, A 2 (s) = A 2 (s)/N 
reduces to the integral 



Ao(s) 



2 r 



7T JO 



dx sin x 



1 — 4p sin 



which for large s behaves as 



A 2 (s 



1 1 



2v^F 



(ps) 



-3/2 



(A8) 



(A9) 



This illustrates the long range of the correlations of the stochastic process n(0, s). 



3- A 3 (j» 

Eq. (3.16) defines 

A *M = ^ + l%<™(m)™(t{ j + l)) ■ < A10 > 

We mainly are interested in the sum 

^ao=p(^+E^o\«)) • ( An ) 

Straightforward summation yields 

I , • n Pt N 1 1 f (, 1 \ • 
A 3 (?, t) = — + + + - 1 j 

dU ' 7 iV 3 2 6iV 2iV V 2iV/ J 

1 ^ , CQS(#) COS + 
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where we used the relation 

£~(£)-(SH)H 

to simplify the result. Comparing to Eq. (A. 2) we find 

A 3 (0,t) = A 1 (0,t) , 

as it should. For iV — > oo, i,j fixed, A 3 (j,t) = -^A 3 (j,t) takes the form 

i /- - 1 i 2 - 2 cosinkj) __.ai.2F / 1 . i\ 

^ 3 (j,t)=t+3 + y-J-^ E — . (A14) 



Thus 



i 3 (j, t) _4i(0, t) + const = AjQ, t) + const . (A15) 



For j > 0, t — > 0, A 3 vanishes exponentially. 



i3(j,t) f ~4^ e " ?/4f • (A16) 



APPENDIX B: SPITZER'S IDENTITY AND END BEAD MOTION 



For an uncorrelated discrete stochastic process n s , Spitzer (17| has derived an expression 



for the generating function of the simultaneous distribution of n t and 

n max {t) = max n s . (Bl) 

se[0,t] 

(It is understood that no = 0, by definition.) For a process symmetric in its hopping 
probabilities, the result reads 



^ ^ ^^max (t) /jjj^max (t) Tit 

t=o 

oo \s 

(1-A)- 1 exp{E — (^ ns +^ ns -2)6(n s -l)} . (B2) 



Taking the derivative with respect to z and putting z = 1 = w, we immediately find 

oo oo \s 

E A*n^(t) = (l-A)- 1 E —n a S(n 8 -l) . (B3) 

t=0 s=l S 

Inverting the discrete Laplace transformation we find the result (5.1): 
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In the same way we can construct an approximation for the fourth moment 

,- v ^ 

gi(Q,N,t) = (roAt)-roAU)) 4 ) ■ (B4) 

Following the argument of Sect. II. A we find the expression 



fi(0,N,t) = (2e s n max (t) -4n(0,t)) 2 

= % (4nl ax (t) - 4n max (t)n(0, t) + n*(0, tj) . (B5) 

Again replacing n(0, s) by an uncorrelated process n s , we use Eq. (B 2) to evaluate this 
expression. We find 



t — 9 t— ] I II I 

ni raJ n,/ 



nl ax (t) = Eg+E (B6) 

=1 s,s'=l 



1 ni 



M)n t = Y.± (B7) 

s=l Z6 



This results in the approximation 

1/2 

IZi l;iM) v; l;iM) «M I _ _ 

g 1 (0,N,t)=£ l 



\n(0,s)\ \n(0,s') 



(B8) 



s,s'=l SS 

Of particular interest is the universal region where all discreteness corrections have died out. 
Eqs. (3.24), (3.26) yield 

n 2 (j,t) = 2p A 1 (j,t) 

2 1/9 

\n(j,s)\ = —j= (poA^s)) 

\/7T 



For times large compared to T2 we can use the approximation (cf. Eq. 5.23) 

Ai(0,s) = Ns , 

and we can evaluate the summations in Eq. (B.8) as integrals. This yields 



- l 3 n » 1.085, T 2 «t«T 3 . (B9) 



<?i(0,A^) V4 
For T < t < T 2 we find (cf. Eq. 5.23) 

Ai(0,s) = 2N ' " 

which in the same way leads to 

^i(0, TV,*) 1 /7T , 2 ^ /INN 1 / 4 



This should be compared to the corresponding result for the central segment 
9i{^,N,t) V 2 

For the end segment the ratio <?i/<?i in all the universal time region is much closer to 1. 
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APPENDIX C: APPROXIMATE FORM OF 

'PmaxJ iP"tni t) 

To construct an approximation for the simultaneous distribution of n max (t) and n(j,t) 
we start from the factorization (5.40) to write 

+00 

]t) = / Pmax,o( n m\n t ;t)V2(n t ,n;0,j,t) (CI) 

nt=— 00 

A rigorous expression for P 2 (^t, 0, j, t) is given in Eq. (3.30). For the conditional probabil- 
ity 'Pmax,o( n m\nt', t) to find n max (t) = n m for given n(0, t) = n t we take the result of a simple 
random walk, making steps ±1 with probability po and steps with probability 1 — 2p Q . 

Let the walk start at no = 0, reach n max (t) = n m at time s for the first time, and end at 
nt < n m . The contributions of the parts [0, s — 1] or [s, t] of this process can be written in 
terms of powers of appropriately truncated matrices of type W (Eq. 3.7), and the evaluation 
is straightforward. We find 

A (rati t)V ma x,o(n m \n t ; t) = 6 (n m - ^(n t + \n t \)^j (1 - S nm0 )p 

4 roo 

— - / dx x dx 2 sin(xi(l + n m + n t )) sin(xi) sin(x 2 ) sm(x 2 n m ) 

7T Z Jo 

^^(xx) - a t ^ l (x 2 ) 

•ao(aJi) 7 — x r^— 

a (xi) - a {x 2 ) 

2 r 00 ,1 
+Q(n m - l)S nmnt po- / dx sm(x)sin(xn m )a (x) 

TC JO 

2 f°° 

+8n m o<$>(-nt)- / sin(x(l - n t )) sm(x)a (x) , (C2) 

7T JO 



where 



ao(^) = l-4j9 sin 2 Q . (C3) 



We now may adjust p so as to reach the desired value of n max {t): 



+00 



n, n p x {n t \t)V maxfi {n m \n t \t) = n max (t) (C4) 

n t =— oo n m =0 

This construction of V max j{n m ,n;t) guarantees that 
i) the relations < n m and n t < obeyed, 



ii) n 2 (j,t) and n max (t) take their proper values, and 

iii) g-i(j,N,t) (Eq. 2.13) evaluated with this distribution for j > 0, t < T R (j) reduces to 
the results of Section IV, and for j = yields gi(0, N, t) as evaluated in Sect. V.B.. 
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Being interested mainly in beads j not extremely close to a chain end (J ^ 20, say), we 
note that pT^{j) ^> 1, so that we may simplify the result by taking the limit of large time: 
olq{x) — > exp(— p tx 2 ). Then the integrals in Eq. (C.2) can be evaluated further to yield 



where we use the notation 



© (n m -\{n t + \n t 



exp 



7T 



4 + 21 



ZmZt 



+ 



(C5) 



n r 



z t 



n ma x(t) y Pot 2 
n t n t 



(C6) 



The neglected terms of relative order yl/p t are due to boundary effects (terms 5 nm0 , S nmTH 
in Eq. (C.2)). 

We now evaluate V-Ant, n; 0, j, t) in the corresponding approximation, which implies ex- 
panding the argument of the exponential in Eq. (3.30) to second order in ip, ip t . The resulting 
Gaussian integral yields 



V 2 (n t ,n-0,j,t) = ^-(n 2 (0,f)n 2 0',0 



2tt 
exp 



-1/2 



2\-l/2 



2 -2a 2 



(z 2 + z 2 — 2a z z t ) 



where 



n 



(n 2 (j,t))V2 ' 



(C7) 



(C8) 



_ i 3 (j,*) 

(^flM)^,*)) 1 ' 2 ' 1 j 

Inserting Eqs. (C.5), (C.9) into Eq. (C.l) and carrying through the summation over n t as 
integral over z t we find the result given in Eq. (5.46). 
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Figure captions 

SORRY: Figures 3 and 5 have an unconventional format and are not included in this 
file. For copies of these figures please contact ebert@lorentz.leidenuniv.nl. 
Fig. 1 

The Evans-Edwards model in 2 dimensions: a piece of a chain on a square lattice, show- 
ing two hairpins. The crosses indicate the obstacles. The hairpins are the mobile units 
corresponding to the 'defects'. 
Fig. 2 

Reptational motion. The full line gives the initial chain configuration. The broken lines are 
the new parts found at time t. Arrows point to j < or j > , respectively. The motion of a bead 
in the original tube is indicated, by its initial (heavy dot) and final (shaded circle) positions. 
Fig. 3 

logi„(0i(j, *)/&«(*)) as Unction of log 10 t for j = 1/100, 1/10, 1/2. g ass (t) = 2£ s Jf (|) 1/4 
represents the first t 1 / 4 -regime. Short dashes: height of the second (a/2-) t 1//4 -regime. Long 
dashes: t 1 / 2 -behavior. 
Fig. 4 

g 1 (N/2,N,t)/g ass (t) as function of \og 10 (pt), where g ass (t) = (t/vr) , cf. Eq. (4.13). 

Full line: N — > oo, long dashes: N = 40, 100, 1000, short dashes: The t 1//4 -asymptote, here 
normalized to height 1 by extracting g ass (t). 
Fig. 5 

Ratio <7i(0, t)/gi(^, t) as function of log 10 t. 
Fig. 6 

gi (0, N, t)/g ass (0, t) asymptote as function of log 10 (pt). g ass {0, t) is given by Eq. (5.28). Full 
line: N — > oo. Long dashes: N = 40, 100, 1000. The short dashes give the i 1//4 -asymptote 
here found at 1 by normalization. 
Fig. 7 

Schematic plot of g±(j, N, i)/<7i(y, N, t) for 1 <^ j <^ N, showing the sequence of crossover 
behavior discussed in the text. Division by <7i(y, N, t) serves to transform power law regimes 
to plateaus indicated by the thin lines. 

Fig. 8 

Simultaneous distribution of n max (t),n(j,t) for values n max {t) / n max {t) — z m — 1,2,3 as 
function of z — n(j, t)/(n 2 (j, t)) 1 / 2 . Full lines: a(j,t) = 0.9; broken lines: a(j,t) = 0.99. 

Fig. 9 

Results for gi(j, N = 10 3 , t) normalized to the initial t 1 / 4 -behavior g ass (t) = ^ i^/ 71 ) ' 
Eq. (4.13), for a) j = 20; b) j = 100. Fat full lines give the full result. Thin lines: . .). 
Dot-dashed: g r (. . .). Dotted: g(0, 10 3 ,t). The arrows indicate T R (j). 
Fig. 10 

Reptation time as function of chain length. Fat line: full model. Thin line: continuous chain 
limit. Long dashes: T 3 ~ iV 3 . Short dashes: T 3 ~ A^ 3 - 4 . 
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